Instabilities and waves in thin films of living fluids 
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We formulate the thin-film hydrodynamics of a suspension of polar self-driven particles and show 
that it is prone to several instabilities through the interplay of activity, polarity and the existence of a 
free surface. Our approach extends, to self-propelling systems, the work of Ben Amar and Cummings 
[Phys Fluids 13 (2001) 1160] on thin-film nematics. Based on our estimates the instabilities should be 
seen in bacterial suspensions and the lamellipodium, and are potentially relevant to the morphology 
of biofilms. We suggest several experimental tests of our theory. 



Active hydrodynamics, the collective behaviour of self- 
driven, orientable particles in a fluid medium, is a topic 
of intense current research [TJ|U[31l31[51[ni|7JIH]- In this 
Letter we study the dynamics of a film of fluid, thin in 
the z direction and spread on a solid surface in the xy 
plane which is also the easy plane for the orientational 
order parameter p of the polar active particles suspended 
in the fluid. For an active system, this polarity implies 
a current v^cp with respect to the fluid, where is a 
characteristic drift velocity and c the concentration of 
active particles. While our formulation is general, we 
study mainly the properties of perturbations about an 
ordered, uniform reference state (p) = ic. 
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FIG. 1: Top view of film with contractile polar filaments 
(solid arrows) with a preference for pointing downhill. Active 
stresses cause fluid flow (dashed arrows) towards the open 
end of a splay perturbation, leading to gradients in the height 
of the fluid film. Modes propagating forward experience a 
height gradient that further torques (dash-dotted arrow) the 
filaments in the direction in which they were already per- 
turbed. 

Our main results are as follows: (i) Active, ordered 
thin films, although dominated by viscosity, not inertia, 
show a wavelike response to external disturbances, as a 
result of the coupled dynamics of free-surface undula- 
tions the active stress field, and the concentration, (ii) 
In large regimes of parameter space this coupling pro- 
duces a novel instability whose growth rate, to leading 
order in small wavevector k = {kx,ky), varies as kykV"^ 
for kx ^ £,ky and as fcj for k^ ^ ^fc^. The crossover 
length ^ = r|Ccro|/io/wQ^ depends on the drift velocity 



of the active particles, the mean film thickness /iq, 
the viscosity /x, an orientation mobility F, the coupling 
C, present in equilibrium systems as well, of particle ori- 
entation to free-surface tilt, and the typical active con- 
tractile or tensile stress uo = Wcq. W is the strength 
of the force-dipole associated with each active particle, 
and Co the mean concentration of active particles, (iii) 
Low motility, i.e., small Vq, promotes the instability and 
high motility suppresses it, a point we will return to at 
the end of this paper. The instability manifests itself in 
moving (convective) and static (absolute) form. We also 
offer a simple physical picture (Fig. [ij of the mechanism 
underlying this unique instability, (iv) We find several 
other possible modes of free-surface instability, which we 
discuss towards the end of the paper. 

Before presenting our model in detail, a few words are 
in order about the systems to which it might apply. Bac- 
terial colonies on a surface frequently take the form of 
organised, rigid biofilms [3], e.g., those of Pseudomonas 
aeruginosa |10j . An understanding of the formation of 
these structures should begin with the thin-film hydro- 
dynamics of active fiuids. The lamellipodium of crawling 
cells is propelled by the ATP-driven force of a continu- 
ally polymerizing and depolymerizing, contractile actin 
network just below the cell membrane [TT] and is thus a 
moving, thin, active fiuid film. 

We therefore consider a fiuid film containing active 
particles with concentration field c, and orientation de- 
scribed by a vector field p = (p_L,Pz)- Since we are 
interested here in perturbations about a macroscopically 
ordered state, we set |p| = 1. We construct equations of 
motion for c{t±, t), p(rj^, t), and the height field h{r±,t), 
i.e., the film thickness, as functions of in-plane position 
rj^ and time t. Our treatment generalizes [12] to the case 
of active systems. 

The kinematic boundary condition [13) h — — u± ■ 
W±h connects h to the hydrodynamic velocity field u = 
{uj_,Uz) evaluated at the free surface. Incompressibility 
V • u = leads to the simplification 

dth + \/j_- [ uj_dz^O. (1) 
Jo 



talso at JNCASR, Bangalore 560 064 India 



We eliminate Uj^ in favour of p, c and h through the 
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Stokes equation in the lubrication approximation [I'd] 
= 0, |Viu|<|9,u|: 



(2) 



where fj, and P are the viscosity and pressure field of the 
fluid, and cr"(r) = VFc(r)p(r)p(r) is the intrinsic stress 
field Umi] of the active particles. We discuss the role of 
gravity later in the paper. 14^ < and W > correspond 
respectively to contractile and tensile activity. 

The active-particle current is along p, and the particles 
cannot leave the film. Thus Pz{z = 0) = 0, and p - N — 
at z = h, where N = {-V ±h,l)/ ^/T+JV^Jif is the 
outward normal to the free surface. So pz ~ p_l ■ Vj^/i 
at the free surface, interpolating linearly, through di- 
rector elasticity [IS], to = at the substrate, i.e., 
Pz — {z/h)p± ■ \/±h. Therefore in a z-averaged de- 
scription pz = {l/2)dxh and dzPz — h~^dxh. Con- 
sider small deviations pJ6] about a state aligned along 
x: = X + Oij, 9 <^ 1. From the foregoing discus- 
sion, the active force density has components Vicrfj. = 
(yo{dy9 + dxc/cQ + h~'^dxh), Viaf^ = aod^O and Viaf^ = 
CTo9^/i/2 to linear order. We now eliminate the pressure 
P in favour of h, p and c, through the z-component of 
^ and stress continuity at the free surface z = h, lead- 
ing to P{x,y,z,t) = Po — j\/^h + (To{h — z)dlh/2, where 
Pq is a reference pressure and 7 the surface tension of 
the fluid [13]. Since we are interested in the dynamics at 
large scales in the xy plane, we treat the active stresses 
in a z— averaged approximation, allowing us to borrow 
the methods of [13, for spreading under gravity, whose 
force density is 2— independent. 

Using this in ^ and integrating twice over z with 
dzU±{h) = and uj^(O) = we get 

\ . 

to linear order in perturbations of h, where = 

ao[{dye + dxc/cQ + h-^dxh)x + djy]. Using ^ in 
linearizing h ^ ho + Sh, c = cq + Sc, the in-plane fouricr- 
components (5/ik(^), Sck{t), 9^(1) obey 
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The four effects of activity on the right-hand side of 
Q are, from left to right within the square bracket, 
(i) curvature-induced flow; (ii) anisotropic osmotic flow; 
and, in parentheses, (iii) splay-induced flow from tilting 
the free surface and (iv) an active anisotropic contribu- 
tion to the effective tension. The final term on the right 
of ([4]) is ordinary surface tension. For contractile ((To < 0) 
stresses, term (iii) destabilises, and the active tension in 
(iv) stabilises, the surface; for tensile {ctq > 0) stresses the 



opposite happens. The behaviour of the active tension 
term is consistent with the idea that contractile stresses 
pull in along the long axis of the particles, giving addi- 
tional elastic resistance to stretching along that axis. 

The dynamics of the polar orientation field p dif- 
fers from that of the conventional nematic director. 
First, symmetry cannot rule out a coupling of the form 
— C / cPxp± ■ Vj_/i which couples p_L to tilts of the free 
surface |T7]. For C > a tilt of the free surface tries to 
make p_L point uphill. 

We can think of two possible microscopic origins for 
such a term. Consider an imposed thickness variation 
tilt V±h of the free surface. One end of each polar par- 
ticle is in general expected to be fatter than the other, 
and the particles will be best accommodated with the 
fat end oriented towards the direction of increasing h. In 
addition, the spatial variation in the separation between 
the free surface and the fixed substrate should give rise 
to an electric field in the plane. Polar particles in general 
have an electric dipole moment, and will thus be oriented 
by this field. Secondly, consider the spontaneous-splay 
term [IB] Hgp = — J cPxCV ■ p in the effective Hamil- 
tonian for p, where C is a phenomenological parameter 
which can depend on the local concentration, say C = 
C{co) + C'{co)Sc+.. . = C + C'Sc. Through the "equilib- 
rium" part of the dynamics of our system, the C and C" 
terms will contribute -TSH/6p = TCS{z)V j_h-TC'VSc 
to the equation of motion for p, F being a kinetic coeffi- 
cient. On dimensional and symmetry grounds we expect 




FIG. 2: An object with head and tail of sizes £ and £' re- 
spectively, and length £, produces a spontaneous curvature 
1/7? ~ {£-£')/£'^. 

C ^ J- jR, where is a typical force scale and i?, the ra- 
dius of curvature associated with the particle shape (see 
Fig. [2])., is the particle thickness divided by the fractional 
length difference between the two edges of the particle. 

Combining the above with the results of [3], the dy- 
namics of the the angle field 9 becomes 



(4) Q^Q ^ -a^vQdx9-idybc+8{z)VCdyh+(\Ay 



(5) 

The first term on the right-hand side of (|5| is advec- 
tion of B with a speed proportional to [T9j . The rest, 
from left to right, are coupling to concentration gradi- 
ents, with C = FC", coupling to the free surface through 
the coefficient C, ffow-orientation coupling as in nemat- 
ics [15], and gradient elasticity of p. The first three 
terms are polar, and the first is a consequence of activity. 
Aij = i {ViUj + V jUi) and Q^ij = ^ {WiUj — WjUi). For 
stable flow-alignment [T5j |A| > 1. 
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Eliminating u in favour of 9 and h through ([s]), and 
averaging over z, ^ becomes, to leading orders in k and 
linear order in the fields, 
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dt9 = + — ky5h^-{D+kl + D^kl 
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- {i(ky ~ ^kxky)Sc, 



(6) 



where D± = D - {X ± l)hl<7o/4:fi, and $ = (A - 
l)/i^(To/4coAi. 

The active contribution to particle currents with re- 
spect to the fluid is VQcp where vq is a drift velocity. 
Thus in the lab frame of reference the continuity equa- 
tion for the concentration, apart from diffusion, is simply 
dfSc = —V - [(u-|-i;op)c]. Linearizing around cq and using 
([5]) and overall incompressibility, 

dtSc — —icQVokyO — ivokxSc + 0(fc^(5ck, k^Shy^). (7) 

The 0(fc2) terms in ^ come from including fluctuations 
in px as well by explicitly including diffusion of particles. 

The instability we are chiefly concerned with here, aris- 
ing from the combination of activity and the tilt coupling 
C, is best seen in the extreme case of immotile but active 
polar particles, i.e. f = but erg 7^ 0- To leading order 
in k the concentration plays no role, and the dynamics is 
determined by (HI) and dsj), yielding an unstable mode in 
the system which grows at a rate kx ky. In detail, for 
time-dependence exp(~iLut), to leading order in k, the 
mode frequency 



c. = ±l±l^S^|^(^)V^|Caofc.|V^|fc,| (8) 

has a real and an imaginary part; the mode displaying 
an instability travels in the -l-a; direction (—a; direction) 
if (JoC > (croC < 0). 

This instability can be understood from a simple phys- 
ical picture, involving the interplay of active stresses (Tq 
and spontaneous splay C. To fix ideas, consider as in Fig. 
[1] a long-wavelength splay perturbation on a set of con- 
tractile filaments (Tq < aligned with polarisation along 
-|-x, with a preference for pointing downhill (C < 0), in 
an initially fiat film. The active stresses then generate 
a fiow in the film that increases the height of the free 
surface ahead of the outward-splayed filaments, and de- 
creases the height ahead of the inward-splayed filaments. 
The mode propagating along +x will then cause the fila- 
ments to sample a height gradient that tilts the filaments 
in the centre of the field further to the right. The same 
argument, mutatis mutandis, goes through for the tensile 
and/or C > cases. 

For motile particles, 7^ 0, let us for the moment 
continue to ignore the concentration field. In the limit 
TC /vq ^ 1, ^ and ([s]) then yield mode frequencies 
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(10) 



We now include the concentration via ([t]) , but continue 
in the simplifying limit FC/wq 3> 1. One mode retains 
the purely relaxational form (|9| , with coefficient modified 
at 0(C)- The propagating mode in (10 1 becomes a pair 
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with speeds of order wo ± 0(C), with relaxation rates ~ 
±{hl<JoCV / ^lVo)[l + 0{C,)]k'^ . The results from ^ to this 
point establish the main claims at the start of this paper. 
In particular, it is readily seen that ([9|, (|10| crossover to 
(g for kx < £,kl with ^ = T\C(7o\hl/vlti. 

Some general features are worth noting: (a) Regardless 
of the sign of a^C, there is always an instability; (b) in- 
creasing Vq weakens the instability, presumably because a 
collection of filaments drifting along x samples alternat- 
ing height gradients along y whose effect cancels; (c) for 
(ToC > the instability moves with a speed ~ wq ± 0{C,), 
combining the effects of the drift speed and the "pres- 
sure" due to C, while for a^C < there is an instabil- 
ity in mode wi, predominantly the height field h, that 
grows without travelling, the remnant of the instability 
discussed for vq = 0. 

It is crucial to note here that activity and polarity con- 
spire to produce the instability, at leading order in gradi- 
ents, even at zero motility. Motility, also a consequence 
of active polarity, is encoded in the active stress only at 
next-to- leading order [T?]. The stabilising effect of in- 
creasing Vq at constant ctq is thus not paradoxical. 

For TC/vq ^ 1 all instabilities involving the interplay 
of activity with the spontaneous-splay coupling C disap- 
pear. The remaining instabilities involve neither polarity 
nor a drift velocity. The most interesting of these, men- 
tioned briefly below arises as follows: Consider a 
flat free surface with contractile filaments aligned along 
X, and impose a small tilt 6h cx x. The filaments at 
the free surface are then tilted relative to those at the 
substrate. The resulting splay in the xz plane implies 
that active stresses will pump fiuid towards the open end 
of this splayed configuration, thus increasing the tilt. In 
addition, the instabilities of a bulk active ordered suspen- 
sion as originally discussed in [31 |S] can arise, modified 
here by confinement to a thin film, and the anisotropic 
active tension [mentioned after (|4|], if large enough, will 
destabilize tensile active films at order k^. 

We now compare the relaxation rates in ([s]) with the 
stabilizing contributions pgh^k^ / ^ due to gravity, p being 
the mass-density of the film and g the acceleration due to 
gravity, ^h^k'^ / p due to surface tension, and Kk^/fi due 
to orientational relaxation, where K is a, Frank constant. 
On dimensional grounds and from standard liquid-crystal 
physics [15] we take F « l//i, where fi « 10"'^ Pa s is the 
viscosity of the medium. With this estimate we find the 
instability survives these stabilizing agencies provided 
kho KmrniCao/p^g^hliCaoho/j^/^Caohl/K^]. For 
7 and p we use values for water, cto should be of order 
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faco = (pf/a^ where / is the force exerted by the active 
particle on the fluid, a the particle size, and we take the 
particle volume fraction = cqq^ of order unity since 
we want an ordered phase. For a bacterium of size a ~ 
a few /im moving at several /itm/s through water / ^ 1 
pN. We argued that C ~ !F/R where T is a. force scale 
and R the spontaneous curvature associated with a po- 
lar object. For thermal hard-rod systems T ~ kBT/£ 
where is a rod length. For /im-sized objects this would 
give C ^ 10~^ dyn/cm. For a self-propelled object, it is 
possibly more reasonable to use the force-dipole strength 
W — af as the natural energy scale, in which case T ^ j . 
For an object of radius 2 /xm moving through water at 
20 jjxn/s this gives J- in the range of a piconewton and 
thus C ^ 10^"* dyn/cm. For film thickness /iq — 20 /j,m, 
and C ~ 10^^ dyn/cm, the instability should be seen 
if fc/iQ ^5 10^'^ . However, this is a pessimistic estimate: 
First, activity rather than thermal energy is quite likely 
responsible for C . Secondly, the tension 7 for biofilms is 
that between the film and the ambient aqueous medium 
and therefore much smaller than the air- water value. We 
expect therefore that the instability should be seen over 
a much wider wavenumber range, possibly fc/io ~ 1- 

The experiment of choice to test our ideas would be 
to compare the dynamics of two initially flat thin fllms 
of uniform concentration, one with highly motile bacte- 
ria, the other with a low- motility mutant, under condi- 
tions of constant bacterial concentration. The mutant 
population would correspond to a system with vq small, 
and should show our instabilities. The observations of 
[lOj on Pseudomonas aeruginosa are of interest in this 
regard, but are complicated by the fact that the bacte- 
ria are dividing. The relation of our mode structure to 
the excitations seen in studies [8] of Myxococcus xanthus 
and to the lamellipodium of crawling cells [TT] remains 
to be explored. A complete treatment must include actin 
treadmilling, as well as the elasticity and anchoring of the 
cell membrane. In all these examples, the dynamics of 
the bulk fluid above the fllm must also be included. Since 
inhomogeneities in orientation give rise to mass flux, the 
instabilities will produce large concentration variations 
which will be important especially when going beyond 
the linearised treatment. 

We thank S. Banerjee for pointing out an error in 
an earlier version of this paper, J. Frost, K. Vijay Ku- 
mar and J. Kas for valuable inputs, and CEFIPRA 
project 3504-2 and the DST Math-Bio Centre grant 
SR/S4/MS:419/07 for support. 



* Electronic address: 'sumithra®phvsics.iisc.er net.in | 

* Electronic address: sriramOphysics. iisc.ernet.ini 

[1] J. Toner and Y. Tu, Phys. Rev. E 58, 4828 (1998). 

[2] J. Toner, Y. Tu and S. Ramaswamy, Ann. Phys. 318, 



170 (2005); F. Jiilicher et al, Phys. Rep. 449, 3 (2007). 
[3] R.A. Simha and S. Ramaswamy, Phys. Rev. Lett 89, 
058101 (2002). 

[4] T.B. Liverpool and M.C. Marchetti, Europhys. Lett. 69, 
846 (2005); Phys. Rev. Lett. 90, 138102 (2003) 

[5] K. Kruse et al, Phys. Rev. Lett. 92, 078101 (2004). 

[6] A. Ahmadi et al, Phys. Rev. E 74, 061913 (2006); L. 
Giomi et al, arXiv:0805.1680; T. B. Liverpool, M. C. 
Marchetti, Phys. Rev. Lett. 97, 268101 (2006); J.-F. 
Joanny et al, New J. Phys. 9, 422 (2007); S. Ramaswamy 
and M. Rao, New J. Phys. 9, 423 (2007); R. Voituriez, 
J-F Joanny and J. Prost, Europhys. Lett. 70, 404 (2005) 
D. Saintillan and M. J. Shelley, Phys. Rev. Lett. 100, 
178103 (2008); V. Mehandia and P. Nott, J. Fluid Mech. 
595, 239 (2008); T. Ishikawa and T. J. Pedley, J. Fluid 
Mech. 588, 399 (2007); C. Wolgemuth, Biophys. J. 95, 
1564 (2008); M. E. Gates et al, Phys. Rev. Lett. 101 
(2008) 068102; L Llopis, L Pagonabarraga, Europhys. 
Lett., 75, 999 (2006); S. Ramachandran, P. B. S. Ku- 
mar, L Pagonabarraga, Eur. Phs. J. E 20, 151 (2006); 
J. P. Hernandez-Ortiz et al, Phys. Rev. Lett. 95, 204501 

(2005) ; P. Underhill et al, Phys. Rev. Lett. 100, 248101 
(2008); G. Gregoire, H. Ghate, and Y. Tu, Phys. Rev. E 
64, 011902 (2001). 

[7] G. Dombrowski et al., Phys. Rev. Lett. 93, 098103 (2004); 
L Tuval, et. al, Proc. Natl. Acad. Sci. 102, 2277 (2005); 
T. J. Pedley and J. O. Kessler, Annu. Rev. Fluid Mech. 
24, 313 (1992); L S. Aranson et al., Phys. Rev. E 75, 
040901 (R) (2007); X.-L. Wu and A. Libchaber, Phys. 
Rev. Lett. 84, 3017 (2000). 
[8] O. Ignoshin, R. Welch, D. Kaiser and G. Oster, Proc. 

Natl. Acad. Sci. 101, 4256 (2004). 
[9] http: / /en. wikipedia.org/ wiki/Biofilnil 
[10] M. Klaussen et. al, Mol. Microbiol. 50, 61 (2003); ibid. 

48, 1511 (2003). 
[11] K. Kruse et al, Phys. Biol. 3 (2006) 130137 
[12] M. Ben Amar and L.J. Gummings, Phys. Fluids. 13 1160 
(2001) 

[13] H.A. Stone, in Nonlinear PDEs in Condensed Matter 
and Reactive Flows, NATO Science Series G: Mathemat- 
ical and Physical Sciences, 569, H. Berestycki and Y. 
Pomeau eds., Kluwer Academic, Dordrecht, The Nether- 
lands (2002); A. Oron, S.H. Davis and S.G. Bankoff, Rev. 
Mod. Phys. 69, 931 (1997). 

[14] Y. Hatwalne, S. Ramaswamy, M. Rao and R.A. Simha, 
Phys. Rev. Lett 92, 118101 (2004). 

[15] P.-G. de Gennes and J. Prost, The Physics of Liquid 
Crystals, Glarendon, Oxford (1993). 

[16] Even in a linearized description, px is not constant but 
is slaved to the concentration and height fields, as seen 
by generalizing the treatment of PJ to include the effects 
of the free surface. This results in subdominant contri- 
butions to the dynamics of the height field Q, which we 
ignore here. 

[17] We thank Shiladitya Banerjee for pointing out the error 
in our arguments leading to this term in an earlier version 
of this paper. 

[18] W. Kung, M.G. Marchetti, K. Saunders Phys. Rev. E 73 

(2006) 031708. 

[19] The proportionality constant a\ is not necessarily unity. 
Since the velocities are defined with respect to the fluid's 
frame of reference, shifting all velocities by a constant is 
not an invariance pPj. 



